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It has been suggested that Glauber (inflow) and Sznajd (outflow) zero-temperature dynamics 
for the one dimensional Ising ferromagnet with the nearest neighbors interactions are equivalent. 
Here we compare both dynamics from analytical and computational points of view. We use the 
method of mapping an Ising spin system onto the dimer RSA model and show that already this 
simple mapping allows to see the differences between inflow and outflow zero-temperature dynamics. 
Then we investigate both dynamics with synchronous, partially synchronous and random sequential 
updating using the Monte Carlo technique and compare both dynamics in terms of the number of 
persistent spins, clusters, mean relaxation time and relaxation time distribution. 

PACS numbers: 05.50.-|-q 



I. INTRODUCTION 



The majority of natural phenomena observed in 
physics, biology, geology, social sciences, etc., are non- 
equilibrium processes. Unfortunately, the theory of 
non-equilibrium statistical mechanics is far less devel- 
oped than its equilibrium part. As a result, the most 
ubiquitous phenomena are poorly understood The 
zero-temperature dynamics of simple models, such as 
Ising ferromagnets, provides interesting examples of non- 
equilibrium dynamical systems with many attractors 
(absorbing configurations, blocked configurations, zero- 
temperature metastable states) 2]. In this paper we fo- 
cus on the so-called single spin flip dynamics for the one 
dimensional Ising ferromagnet. The best known exam- 
ple of such a dynamics for the Ising model is Glauber 
dynamics It can be viewed as "inflow" dynamics, 
since the center spin is influenced by its nearest neigh- 
bors j^. Another type of dynamics, which can be called 
"outflow" dynamics, since the information flows from the 
center spin (or spins) to the neighborhood, has been in- 
troduced to describe opinion formation in social systems 
. It has been suggested 0, Q that both dynamics for 
the Ising ferromagnet with the nearest neighbors interac- 
tions are equivalent, at least in one dimension. However, 
it seems to be true only in some particular cases. The 
aim of this paper is to compare generalized outflow and 
inflow dynamics for a chain of Ising spins and show in 
which cases there are equivalent and in which they dif- 
fer. 

In the first two sections we recall ideas of inflow and 
outflow dynamics and formulate the generalized versions 
of both dynamics. We take both dynamics under a com- 
mon roof reformulating them without using the concept 
of energy. In the third section we use the illuminating 



method of mapping an Ising spin system onto the dimer 
RSA model and make simple mean field like calculations 
to show the difference between the two dynamics. In the 
fourth section we present Monte Carlo results for both 
dynamics with several kinds of updating, including syn- 
chronous, partially synchronous and random sequential 
updating. The summary of simulation results and con- 
clusions are the subjects of the fifth and sixth section of 
the paper. 



II. INFLOW DYNAMICS 

The best known example of such a dynamics for the 
Ising model is Glauber dynamics. Within Glauber dy- 
namics, in a broad sense, each spin is flipped Si{T) 
—Si{T + 1) with a rate W{6E) per unit time and this 
rate is assumed to depend only on the energy difference 
implied by the flip 2]. The two most common choices 
of flipping rates in the case of discrete updates are heat- 
bath and Metropolis, both obey the detailed balance con- 
dition: 



W{SE) 



(1) 



Recently it was shown that there was a vast fam- 
ily of dynamical rates, besides these two choices, which 
obeyed condition (Q. Among them the class of zero- 
temperature dynamics defined as: 



W{SE) 



1 USE < 0, 
Wo ifSE^O, 
ii6E> 0. 



(2) 
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occurred to be very interesting. The zero-temperature 
limits of the heat-bath and Metropolis rates are W^^ = 
1/2 and = 1, respectively. For any non-zero value 
of the rate Wo corresponding to free spins, the dynamics 
belongs to the universality class of the zero-temperature 
Glauber model. This is a prototypical example of phase 



2 



ordering by domain growth (coarsening) . The typical size 
of ordered domains of consecutive | and | spins grows as 
L{t) ^ t^/^. The particular value Wq = corresponds 
to the constrained zero-temperature Glauber dynamics 
(0 and references therein). In the constrained zero- 
temperature Glauber dynamics, the only possible moves 
are flips of isolated spins and the system therefore even- 
tually reaches a blocked configuration, where there is no 
isolated spin 0- Very interesting results for the zero- 
temperature Glauber dynamics have been also obtained 
using computer simulations |3. 1^ ITol . 

In out of equilibrium systems, there is usually no en- 
ergy function and the system is only defined by its dy- 
namical rules 12]. This is also the case of the socio- 
physics Sznajd model. For this reason we reformulate 
the definition of the zero-temperature Glauber dynamics 
for the Ising ferromagnet, without using the concept of 
the energy, in the following way: 



1 

— S'i(r) with prob 
-1 



Wo 



if Er. 

if 

if Er 



>0, 
= 0, 
< 0, 



(3) 



where Xlnn ^nn denotes the sum over nearest neighbors. 

In one dimension, which is the case of this paper, the 
above definition can be written as: 



S^{T+l) = 

1 if ^(t),_i + ^(t),+ i > 0, 

-S^{t) with prob Wo ifS'(T)i_i 
-1 if ^(t),-i 



S{t),+i < 0. 



If Si{T)Si+i{T) = 1 then Si-i{T -f 1) ^ S^{T) and 

Si+2{T -I- 1) ^ S^+i{t) 

is equivalent to the rule: 

if S,{t)S,+i{t) = 1 then S,_i(t + 1) ^ S,+i{t) i 
5,+2(T + l)->5.(r); 

while the antiferromagnetic rule states that 

if S,{T)S,+iiT) - -1 then S,-i{t+1) ^ S,+i{t) and 

Thus, the two above rules can be rewritten as a simple 
single rule: 5*1-1 (r + 1) ^ '5'.i+i(t) and S'.j+2(r + 1) ^ 

In later papers [l^, Hyl we have proposed two modifi- 
cations of the model in which the antiferromagnetic rule 
was replaced by one of rules described below: 

modification 1: If S',(r)S'i+i (t) = -1: 5'j_i(t+ 1) ^ 
Si-i{t) and S'i+2(r + 1) £'1+2(1'). 



or 



modification 2: If S',(r)S'i+i (t) = -1: S^-i{t + 1) 

-S'i_i(r) and Si+2iT + 1) -Si+2{T) with prob- 
ability 1/2. 

A generalized dynamics, which includes the two above 
modifications, can be written as: 



^.(r + l) = 
1 



if S,+i(r) + 5,+2(r) >0, 



-S^ir) with prob Wo if (r) + 5,+2(r) - 0, 
-1 if 5,+i(r) + 5,+2(r) <0. 



(5) 



III. OUTFLOW DYNAMICS 



It is easy to notice that modification 1 corresponds to 
= and modification 2 to ^ 1/2. 



The outflow dynamics was introduced to describe opin- 
ion changes in the society. The idea was based on the fun- 
damental social phenomenon called " social validation" . 
However, in this paper we do not focus on social appli- 
cations of the model (interesting reviews can be found in 
CI El El El)- On the contrary, here we investigate the 
dynamics from the theoretical point of view. 

In the original model a pair of neighboring spins Si and 
Si+i is chosen and if SiSi+i — 1 the the two neighbors 
of the pair follow its direction, i.e. Si-i — + Si{= S'i+i) 
and iS'i-i-2 Si{— Si+i). Such a rule has been used also 
in all later papers dealing with the one dimensional case 
of the model. However, the case in which SiSi+i = —1 
is far less obvious. For example, in the original paper 
SiSi+i = —1 led to Si-i — > S'i+i and S'i+2 ^ Si. How- 
ever, some authors have argued lU] that such a rule is 
unrealistic in a model trying to represent the behavior of 
a community. Moreover, the original Sznajd model with 
both the ferromagnetic and the antiferomagnetic rules is 
equivalent to a single simple rule that every spin takes 
the direction of its next nearest neighbor independently 
of the product SiSi+i. To see this observe that the fer- 
romagnetic rule: 



IV. MAPPING ONTO THE DIMER MODEL 

As was mentioned in the previous sections, for Wq — 
the system under inflow (Glauber) dynamics described 
by formula l@J eventually reaches a blocked configura- 
tion, where there is no isolated spin. On the other hand, 
the system under the outflow dynamics described by (O 
always reaches the ferromagnetic steady state. Thus, for 
14^0 — the difference between the outflow and inflow 
dynamics is obvious. Nevertheless, within the mean field 
approach |23| and the Galam's unifying frame 6] both 
dynamics are equivalent, i.e. there is no difference be- 
tween the outflow and inflow dynamics, even for Wq = 

Here we use the illuminating method of mapping the 
Ising spin system onto the dimer RSA model. This has 
been already done for the inflow dynamics 0: 

j Xi — SiSi + l = 1 => O, 

\X, = SA+i = -1 ^ .. 

In the case of inflow dynamics the following transitions, 
which change the state of the system, are possible: 
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spins 


particles 


J- T i- ^ J- J- J- 
TiT^TTT 


• • — > oo 

• • — > oo 


Itt^l It 

•4' 1 1 -1- -i- 1 


Wo 

• O — > 09 

Wo 

• o — > o» 


TTi-tiJ, 


Wo 
om — > #0 

Wo 

0* — > 90 



Thus, after the mapping there are only two types of tran- 
sitions for the inflow dynamics: •• — > oo and o» <5 •o. 
This mapping shows at once that for Wq = the dynam- 
ics is fully irreversible, in the sense that each spin flips 
at most once during the whole history of the sample. 

It should be noticed that if we map the system under 
the outflow dynamics onto the dimer model we have to 
take into account four particles, because changing the 
border spin influences the next particle. In this case four 
types of transitions are possible: o»o — > oo», o»» — > ooo, 

• o • S • • o and • • • S #00 (to make it clearer the 
flipped spins are denoted by double arrows in the table 
below) : 



spins 


particles 


rm^tm 


o • o — > o o • 
o • o — > o o • 




o • • — > o o o 
o • • — > o o o 




Wo 

• o • — > • • o 

Wo 

• o • — > • • o 




Wo 

• • • — > 9 

Wo 

• • • — > 9 




Wo 

• • o — > • o • 

Wo 

• • o — > • o • 




Wo 

moo — > • • • 

Wo 

• o o — > • • • 



This mapping shows that for Wq = the outflow dy- 
namics consists of two processes - diffusion of • particles 
in the sea of o empty sites and annihilation of •• pairs. 
Thus our model for Wq with random sequential updat- 
ing reduces to an analytically solvable reaction-diffusion 
system A + A ^ (denote by o the empty place and • 
by the A particle). 

For Wq > we can also use the mean field approach 
(MFA). The mean field results for the outflow dynamics 
without dimer mapping can be found in |20j |. Within 
dimer mapping we take into account correlations between 
pairs of the nearest neighbors. Thus if we apply MFA to 
the mapped system we can expect more correct results 
than obtained and within MFA without mapping. 

Let us denote the number of • particles by Nh and 
define b = In our case, in one time step r, only two 
events are possible - the number of • particles decreases 
hy 2/N with probability 7(6) or remains constant. 



For the inflow dynamics: 

7™(6) ^b^, (7) 

and for the outflow dynamics: 

7°"*(6) ^ {1 ~ b)b^ + Wob^ ^ b^l - b{l ^ Wo)] . (8) 

It is seen that the above results are not precise, since 
there is no dependence between 7™ (6) and Wq for the 
inflow dynamics and the only stable steady state in this 
case is 6 = 0, i.e. the ferromagnetic state, which is true 
as long as Wq > 0. However, as it has been noticed 
this result is not correct for Wo = 0. The same results 
can be obtained using the mean held approach without 
mapping. 

However, for the outflow dynamics MFA with dimer 
mapping give s better results than the basic MFA pre- 
sented in [231 • For Wq = there are two steady states: 
6 = 0, i.e. the ferromagnetic state and 6=1, i.e. the 
antiferromagnetic state. For 6 7^ and 6^1, 7°"* (6) > 
which implies that 6 = is an unstable steady state, 
while 6 = 1 is a stable steady state. This result is in 
agreement with computer simulations l^. For Wq = 1 
there is only one ferromagnetic steady statejwhich was 
also confirmed by the computer simulations 

The differences between the outflow and inflow dy- 
namics can be already seen if we apply the mean fleld 
approach with mapping of the pairs of spins into single 
particles. In the next section we present simulation re- 
sults which show yet more differences between these two 
dynamics. 



V. SIMULATION RESULTS 

The spin updating within both dynamics can be se- 
quential or parallel. In this paper we compare both dy- 
namics for random sequential updating, parallel updat- 
ing and partially parallel updating. From now on we 
call the latter case c-parallel updating. Within this up- 
dating a randomly chosen fraction c of spins is updated 
synchronously. Of course, c = 1 corresponds to parallel 
updating and c = to random sequential updating. 



A. Number of persistent spins 

One of the main quantity of interest in the non- 
equilibrium dynamics of spin systems at zero tempera- 
ture is the fraction of spins, P(t), that persist in the 
same state up to some later time t = Nt (measured in 
Monte Carlo Steps) 21, 22J. In this paper we measure 
the fraction of persistent spins for both the outflow and 
inflow dynamics with c-parallel updating for different val- 
ues of c. The initial conflguration consists of a randomly 
distributed fraction p+(0) of up spins. The number of 
persistent spins for the outflow dynamics with Wq = 
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and random sequential updating has been already investi- 
gated by Stauffer and Oliveira (23| who found agreement 
with the results for the inflow dynamics, i.e. decays with 
time t as l/t~^^^. However, it was found that in higher 
dimensions the exponents for the inflow and outflow dy- 
namics are different [2^. Here we investigate the case 
of the Ising spin chain more carefully, i.e. for different 
values of Wq and c. 

The first difference between the inflow and outflow dy- 
namics can be already seen for random sequential up- 
dating, i.e. c = 0. For both dynamics the number of 
persistent spins decays initially as a power-law ~ t^^ . 
However, for the inflow dynamics the exponent is inde- 
pendent of Wq as long as Wo > 0, while for the outflow 
dynamics the exponent is Wo-dependent 9 = 9{c) (see 
Fig.l). Moreover, for the inflow dynamics the power-law 
describes properly the decay of the number of persistent 
spins for all times, while within the outflow dynamics it 
is valid only for t smaller than a certain value of time 
t*iWo), dependent on the flipping probability Wq. For 
Wq — > we obtain t*(Wo) — > 00 and the evolution of the 
number of persistent spins is the same for the outflow and 
inflow dynamics which is in agreement with the results 
obtained by Stauffer and Oliveira p^ . 

More differences can be seen for partially synchronous 
updating with c > 0. At each elementary time step r 
a fraction c of spins is chosen randomly and is changed 
synchronously. In such a case we have noticed that the 
number of persistent spins still decays with a power law 
for the inflow dynamics. However, for the outflow dy- 
namics the power law is no longer valid. The number of 
persistent spins decays very fast in this case (see Fig. 2). 

We may conclude this subsection with the following 
- the number of persistent spins is c sensible only for 
the outflow dynamics. For Wq > and any value of c 
the number of persistent spins in the inflow dynamics is 
described by a power law with nearly the same exponent. 



B. Number of clusters 

Probably the most natural way to investigate the relax- 
ation process of the consensus dynamics is to look at the 
number of clusters in time. A cluster consists of a group 
of spins, each of which is a nearest neighbor to at least one 
other spin in the cluster, with all spins having the same 
orientation. With such a definition consensus is reached 
when only one cluster is present in the system. For both 
infiow and outfiow dynamics with c-parallel updating the 
number of clusters monotonically decays as t^^/^ for any 
value of c. This result shows that the number of clusters 
in time, although a very intuitive and natural measure of 
relaxation, is not a good quantity for comparison of the 
dynamics. 




FIG. 1: The change in time of the number of persistent spins 
on a chain of A*' = 300 spins for random sequential updat- 
ing (i.e. c = 0) for the inflow (top) and outflow dynamics 
(bottom). 



C. Mean relaxation time 

The differences between the dynamics can be observed 
clearly if we look at the mean relaxation time as a func- 
tion of the initial fraction of randomly distributed up- 
spins p+(0). Within 0-parallel updating (i.e. random 
sequential updating) the relaxation is much slower for 
the inflow dynamics then for the outflow dynamics. This 
is also true for the c-parallel updating with small c. On 
the contrary, within 1-parallel updating (i.e. synchronous 
updating) the relaxation the under outflow dynamics is 
slower then under the inflow (see Fig. 3). 

In general, the relaxation times decay with Wq growth, 
but the dependence between the mean relaxation time 
and Wq is different for the outflow and inflow dynamics. 
Two examples for c = 0.2 and c = 0.5 for several values 
of Wq are shown in Figs. 4 and 5, respectively. It can 
be noted that, for example, for c = 0.5 and Wq — 0.8 
the dependence between the mean relaxation time and 
the initial concentration of up-spins p+(0) is nearly the 
same. 

In Figure 6 we have presented the dependence between 
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FIG. 2: The change in time of the number of persistent spins 
on a chain of A'' = 300 spins for partially synchronous updat- 
ing with Wo = 1/2. The top panel presents results for c = 
and the bottom for c — 0.2. Obviously for c > the number 
of persistent spins decays very fast and cannot be described 
by a power law. 



FIG. 3: The mean relaxation times from the random initial 
state consisting of p+(0) randomly distributed up spins for 
Wo — 0.2. The top panel corresponds to synchronous updat- 
ing c = 1, and the bottom panel to c = 0.2. It can be seen 
that the relaxation under the outflow dynamics is slower then 
under the inflow for synchronous updating. On the contrary, 
the relaxation is much slower under the inflow dynamics then 
under the outflow dynamics for small c. 



the mean relaxation times from the random initial state 
consisting of 50% randomly distributed up spins (maxi- 
mal waiting time) and the flipping probability Wq for the 
inflow and outflow dynamics. It is seen that the depen- 
dence on c is much stronger for the outflow dynamics. 
For the inflow dynamics the mean relaxation time is al- 
most the same for all values of c. On the other hand for 
a given value of c the dependence between < r > and Wq 
is stronger for the inflow dynamics. 



D. Distribution of waiting times 

In paper the mean field approach for the outflow 
dynamics with Wq = has been presented and the dis- 
tribution of waiting times needed to reach the stationary 
state has been found. Recall that for the (5-initial condi- 
tions the distribution of waiting times has an exponential 



tafl 113: 

P,P(T)«^(l-m2)e-2-,r^oo. (9) 

Monte Carlo simulations confirm this prediction both on 
a complete graph and on a chain. In this paper we have 
checked also the distribution of waiting times for different 
values of Wq and c both for the outflow and inflow dy- 
namics. It occurs that the distribution of waiting times 
has an exponential tail for any value of Wq and c, al- 
though the exponent depends on these parameters. An 
example for c = 0, showing the comparison between in- 
flow and outflow dynamics, is shown in Fig. 7. 

VI. CONCLUSIONS 

It has been suggested 0, 0] that the zero-temperature 
outflow and inflow dynamics for the Ising ferromagnet 
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FIG. 4: The mean relaxation times from the random initial 
state consisting of p+(0) randomly distributed up spins for 
c — 0.2 for the inflow (top) and outflow (bottom) dynamics. 



FIG. 5: The mean relaxation times from the random initial 
state consisting of p+(0) randomly distributed up spins for 
c = 0.5 for the inflow (top) and outflow (bottom) dynamics. 



with the nearest neighbors interactions are equivalent in 
one dimension. However, this is certainly not true for 
Wq = 0. This particular value corresponds to the con- 
strained zero-temperature Glauber dynamics where the 
only possible moves are flips of isolated spins and the 
system, therefore, eventually reaches a blocked configura- 
tion, where there is no isolated spin |2| . This can be also 
easily shown using the method of mapping the Ising spin 
system onto the RSA dimer model. On the other hand, 
the outflow dynamics leads to the ferromagnetic steady 
state for any value of Wq. This observation motivated 
us to compare both dynamics more carefully. We have 
performed Monte Carlo simulations for both dynamics 
using random sequential updating, parallel updating and 
c-parallel updating (a randomly chosen fraction c of spins 
is updated synchronously). We have measured, for dif- 
ferent values of Wq and c, the distribution of waiting 
times, the mean waiting time, the decay of the number 
of clusters and persistent spins in time. 

It occurs that the qualitative difference between the 
inflow and outflow dynamics is not visible neither in the 
number of clusters in time nor in the distribution of wait- 
ing times. However, it should be noticed that the relax- 



ation time is different for both dynamics. Nevertheless, 
for both dynamics the distribution of waiting times has 
an exponential tail and the number of clusters decays as 
t^^/^ for any value of Wq > and c. 

Differences between the dynamics appear if we look at 
the dependence between the mean relaxation time and 
the initial concentration of randomly distributed up spins 
for different values of Wq and c. For c — 0, which cor- 
responds to random sequential updating, the mean re- 
laxation time is shorter for the outflow dynamics (e.g. 
for Wo — 0.2 and po — 0.5 it is about 10 times shorter) 
than for the inflow. The mean relaxation time < r > 
decreases with Wq increase for both dynamics, but the 
dependence between < r > and Wq is different for the 
outflow and inflow dynamics. Generally, the mean relax- 
ation time decays faster with increase Wq for the inflow 
dynamics for any value of c. Moreover, with increase c 
the dependence between the mean relaxation time for the 
inflow dynamics and the outflow dynamics vanishes. As 
the results for some values of c and Wq (e.g. c = 0.5 
and Wq = 0.8) the dependence between the mean relax- 
ation times and the initial concentration of up spins is 
identical. Of course, this suggests that for some values 
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FIG. 6: The dependence between the mean relaxation times 
from the random initial state consisting of 50% randomly dis- 
tributed up spins and the flipping probability Wo for the in- 
flow (top) and outflow (bottom) dynamics for different values 
of c. It is seen that the dependence on c is much stronger for 
the outflow dynamics. On the other hand, for a given value 
of c the dependence between < t > and Wo is stronger for 
the inflow dynamics. 



of random sequential updating the relaxation under the 
outflow dynamics is much faster than under the inflow 
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FIG. 7: Probabilities of reaching the steady state in time 
larger than r on a chain of L = 200 spins. The distribution 
of waiting times has an exponential tail for both dynamics - 
inflow (top) and outflow (bottom). 



of parameters Wq and c the relaxation under outflow dy- 
namics is faster than under the inflow dynamics. This is 
indeed true. For c = 1 (parallel updating), the relaxation 
is faster under the inflow dynamics for any value of Wq . 

The second quantity which occurs to behave differently 
for both dynamics is the number of persistent spins in 
time. Main differences are seen for partially synchronous 
updating with c > 0. In such a case we have noticed that 
the number of persistent spins decays with a power law 
for the inflow dynamics (like for c = 0). However, for the 
outflow dynamics the power law is no longer valid. The 
number of persistent spins decays very fast in this case. 

Concluding, the inflow and outflow dynamics differ 
very clearly even in one dimension. There is obvious, 
very strong difference for Wq — 0, but also for Wq > 
both dynamics are qualitatively different. In the case 



dynamics. On the contrary, in the case of parallel up- 
dating the relaxation for the outflow dynamics is much 
slower than for the inflow. It should be mentioned here 
that the outflow dynamics with Wq = and synchronous 
updating has been investigated earlier and it has been 
found that in such a case the possibility of reaching a 
consensus is reduced quite dramatically 24]. Also the 
number of persistent spins is c sensible only for the out- 
flow dynamics. For Wq > and any value of c the number 
of persistent spins in the inflow dynamics is described by 
a power law with nearly the same exponent. Generally, 
it occurs that the outflow dynamics is much more influ- 
enced by the type of updating than the inflow dynamics. 
We believe that this result is very important in the vari- 
ous interdisciplinary applications of the zero-temperature 
single-spin flip dynamics. 
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